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We consider a cavity optomechanical system consisting of a Bose-Einstein condensate (BEC) 
interacting with two counterpropagating traveling-wave modes in an optical ring cavity. In contrast 
to the more familiar case where the condensate is driven by the standing-wave field of a high-Q 
Fabry-Perot cavity we find that both symmetric and antisymmetric collective density side modes 
of the BEC are mechanically excited by the light field. In the semiclassical, mean-field limit where 
the light field and the zero-momentum mode of the condensate are treated classically the system 
is found to exhibit a rich multistable behavior, including the appearance of isolated branches of 
solutions (isolas) . We also present examples of the dynamics of the system as input parameters such 
as the frequency of the driving lasers are varied. 
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I. INTRODUCTION 

The optomechanical cooling of mechanical oscillators 
has witnessed considerable progress in the last few years, 
leading to the expectation that a large class of such 
oscillators will soon be cooled to near their quantum- 
mechanical ground state of vibration. In addition to the 
top-down approach where cooling proceeds by mount- 
ing macroscopic oscillators as moving mirrors in an op- 
tical resonator [HQ — often but not always the end- 
mirror of a Fabry-Perot interferometer, there has also 
been increased interest in considering bottom-up situa- 
tions. In that case the mechanical oscillator consists of 
a momentum side mode of an ultracold atomic system 
trapped inside a high-Q optical cavity with fixed mirrors. 
The trapped atoms can be either a thermal sample Q, 
a quantum-degenerate Bose-Einstein condensate (BEC) 
HLI1], or even a quantum-degenerate gas of fermions @- 
In that bottom-up situation the mechanical oscillator(s) 
consist of collective momentum modes of the trapped gas, 
excited via photon recoil 04H|- Specifically, in the case 
of a condensate the intracavity standing-wave field cou- 
ples the macroscopically occupied zero-momentum com- 
ponent of the BEC to a symmetric superposition of the 
states with center-of-mass momentum ±2hk via virtual 
electric dipole transitions. 

It is known that when considering the mechanical ef- 
fects of light on atoms by quantized light fields, there 
are situations where a standing wave does not lead to 
the same diffraction pattern as a superposition of two 
counterpropagating running waves of equal frequencies 
p^4T^ |. This is because in contrast to a standing wave, 
running waves in principle permit one to extract "which 
way" information about the matter-wave diffraction pro- 
cess. This begs the question whether a description of 
the intracavity light field in terms of a standing wave 
is always appropriate to describe the optomechanical ef- 
fects of feeble light fields on ultracold atoms. As a 
first step toward answering this question, this paper ad- 



dresses the somewhat simpler question of understanding 
the difference between classical standing wave and coun- 
terpropagating light fields, that is, the difference in op- 
tomechanical properties of condensates trapped in, say, 
a Fabry-Perot and a ring cavity. One main consequence 
of the presence of two counterpropagating running waves 
is that in addition to a symmetric "cosine" momentum 
side mode, it is now possible to excite an out-of-phase 
"sine" mode as well. In the optomechanics analogy, this 
indicates that two coupled "condensate mirrors" of equal 
oscillation frequencies but in general different masses are 
driven by the intracavity field. This can result in com- 
plex multistable behaviors that we analyze in some detail 
in the following sections. 

The dynamical interaction between BECs and coun- 
terpropagating light fields has a long history. The coop- 
erative scattering of light and atoms in ultracold atomic 
systems, including experimental [l5rll8| and theoretical 
[13-[23[ studies of superradiance in BECs and coherent 
atom recoil lasing (CARL) |24l- [3lJ has been the subject 
of a number of studies, see Ref. [32| for a brief review. 
The work most closely related to our analysis is Ref. [33[ , 
which considers likewise the collective dynamics of atoms 
in a ring resonator, but in the somewhat simplified sit- 
uation where the amplitude of one of the counterpropa- 
gating fields inside the resonator is fixed. 

The remainder of the paper is organized as follows. 
Section II describes our model of a BEC interacting with 
two counterpropagating fields in a ring cavity. It intro- 
duces a description of the condensate that clearly illus- 
trates an optomechanical analogy involving two mirrors 
driven by the intracavity field. In the semiclassical limit 
where the light fields are treated classically, the dynam- 
ics of these two "mirrors" can be understood in terms of 
an interference effect between the propagating fields. We 
further discuss conditions under which one of the mirrors 
can be trapped into a dark state. Section III discusses the 
steady-state properties of the system, and comments on 
the appearance of isolated branches of solutions, or "iso- 
las" [lil, H|[ • We then turn in Section IV to the dynamics 
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of the system, considering in particular the response of 
the condensate to a sweep of the frequency of the exter- 
nal pump lasers. Finally, Section V is a summary and 
outlook. 



II. MODEL 

We consider a BEC of N two-level bosonic atoms with 
transition frequency cu a trapped inside one arm of a high- 
er optical ring cavity supporting two counterpropagating 
traveling-wave fields of frequency oj and driven by two 
coherent pump fields of identical frequency w p and com- 



plex amplitudes rjx and 7/2. The light fields are taken to 
be far detuned from the atomic transition so that the ex- 
cited electronic state of the atoms can be adiabatically 
eliminated. We further assume that the atom trap is suf- 
ficiently soft that the condensate can be assumed to be 
homogeneous along the axis of the resonator. For sim- 
plicity we neglect transverse effects and describe both the 
condensate and the light fields as plane waves, thereby 
reducing the description of the problem to one dimension. 
Neglecting atom-atom collisions and in a frame rotating 
at the pump frequency w p the model Hamiltonian for this 
system is then 
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Figure 1: Excitation of symmetric X c and antisymmetric X s 
collective density side modes. The excitation of these modes 
is due to the interaction between a classically treated zero- 
momentum mode Co — > ■s/No, and the superposition of two 
travelling wave modes ai and &2- 



where ip( x ) an d &i are the bosonic annihilation operators 
for the atomic field and the counterpropagating cavity 
modes i = 1,2, respectively, k — u>/c, A = uj p — uj 
is the pump-cavity detuning, m is the atomic mass, 
Uq = gg/Aa is the familiar off-resonant atom-photon in- 
teraction strength, with g the vacuum Rabi frequency 
and A Q = lo p — uj a the atom-pump detuning. 

The BEC is initially prepared in a macroscopically oc- 
cupied zero-momentum state from which the atoms are 
then scattered into higher momentum side modes ±2£hk 
by the intracavity optical field, where I is an integer. 
For moderate fields the condensate dynamics can be re- 
stricted to the zero-momentum mode and the first two 
side modes, I = 1 0, 0, and the atomic field operator 
can be expanded simply in terms of these modes. Instead 
of a plane- wave expansion, we find it convenient to use a 



sine and cosine basis, so that 
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i>(x) = —j=cq + W-c c cos(2fcx) + \ — c s sin(2fcir), (2) 
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where L is the cavity length and Co, c c and c s are the 
bosonic annihilation operators of the zero-momentum 
mode, the cosine mode, and the sine mode respectively. 
They satisfy the usual bosonic commutation relations 
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With this expansion of the Schrbdinger 



field the Hamiltonian ([T]) reduces to 



H = n(-A + UaN)^2ala. l -J2ih(i 1 *a l -f] l al) 
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+ —j^-{a 2 a x - a|a 2 )(c^c + c ( ! ) c s ), 



(3) 



where ui r = hk 2 /2m is the atomic recoil frequency. This 
Hamiltonian can be further simplified since for weak opti- 
cal fields and large N the depletion of the initial conden- 
sate remains weak. We then treat the zero-momentum 
mode classically, Co — > VN . We further introduce the 
'position' operators 

Xj = (c] +0^/^/2 (4) 

and the corresponding 'momentum' operators 

Pj = - cj)/V2, (5) 
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where j = {c,s}, resulting in the Hamiltonian ([3]) being 
mapped to the optomechanical Hamiltonian 

H = n(-^ + U Q N)^a\ai-Y J ^*ai-'n i a\) 

+ 2hu r (X 2 c + PI + X 2 + P 2 ) 

+ hUoVNX^alcn + a\a 2 ) 

+ ihU VNX s {a\a-L - a\a 2 ). (6) 

In this form, the Hamiltonian provides a simple physical 
picture: the sine and cosine momentum side modes of 
the condensate behave as a pair of mirrors driven by the 
interference between the two counter-propagating intra- 
cavity light fields. As such, they can be thought of as 
a "bottom-up" realization of coupled mirrors driven by 
the radiation pressure. We return to this point shortly, 
but first derive the coupled equations of motion of the 
light-condensate mirrors system. 

In practice the sine and cosine momentum modes are 
coupled to other momentum side modes as a result of 
the presence of a trapping potential that is otherwise ig- 
nored in our discussion. This results in a damping of 
the population of these modes, and associated noise op- 
erators. The resulting quantum Langevin equations ob- 
tained from the Hamiltonian ((6]) are therefore 

X c = 4uj r P c - -yX c + f xc , 

P c = -Auj r X c - Uo\fN{a\ai + a\a 2 ) - jP c + fpc, 
X s = 4uj r P s - jX s + f X3 , 

P s = -4oj r X s - HJqVn ( a\a\ - a\ a 2 ) - jP s + f pa 

161 = U Q VNa 2 (X c - iX s ) - (A + iK)ai + irji + V2ko}^, 

16 2 = UoV~Nai(X c + iX s ) — (A + in)a 2 + ir]2 + \/2~Ka 2 n , 

(7) 

where k is the line width of the ring cavity, A = A — UoN 
is the Stark-shifted cavity-pump detuning, and the side- 
mode damping rate 7 affects both position and momen- 
tum, as discussed in Ref. [9J. All noise operators, f xc , 
a 1 ", etc. are assumed to have zero mean. 

In the remainder of this paper we consider the sim- 
ple limit where all fields, optical and matter-wave, are 
treated classically. A full quantum description of the 
problem will be considered in future work. Taking the 
expectation values of Eqs. (O and factorizing all opera- 
tor products, for example (c^Ac) — > (S2) (X c ), {a\ai) — > 
(al)(ai), etc, yields the classical equations of motion 

X c = -2 1 X c -{l&u 2 r + 1 2 )X c 

-&u r U Q \/N {a^ati +a\a 2 ), (8) 
X s = -2 7 A S -(16^ + 7 2 )X S 

— 4 i ui r UqVn («2 Oil — o>la 2 ), (9) 
iai = U VNa-z(X c - iX s ) - (A + in)ax + ir]i,(10) 
ia 2 = U VNa-i{X c + iX s )- (A + irc)a 2 + w?2,(H) 



where we have set (di) — > cti, (X c ) —> X c , etc. 

Equations (|8]) and ([9]) are particularly instructive in 
that they show that the sine and cosine side modes 
are driven by out-of-phase components of the interfer- 
ence pattern produced by the counterpropagating optical 
fields. In particular, we observe that if these two fields 
have equal phase and amplitude, then the sine mode is 
not excited: it becomes a dark state as a result of the de- 
structive interference between the coherent recoil effects 
from the two counterpropagating fields. Similarly, the co- 
sine mode can become a dark state for out-of-phase and 
equal amplitude intracavity fields, see Fig. 1. We show 
in the next section that these dark states are useful in 
understanding the steady-state properties of the system. 
Interestingly, in high-Q Fabry-Perot resonators it is usu- 
ally appropriate to describe the intracavity light field as 
a standing-wave, a cosine-wave if the origin is chosen as 
the center of the resonators. In that case, only the cosine 
momentum side mode of the condensate is excited, and 
the problem reduces to the familiar BEC optomechanical 
situation of Ref. Q • 

We remark that even in Fabry-Perot interferometers it 
is strictly speaking never completely appropriate to ig- 
nore the other mode, the sine mode of the optical field, 
unless the mirrors are perfectly reflecting (and infinitely 
heavy). This is especially so in the limit of feeble quan- 
tized fields, certainly the quantum fluctuations of the co- 
sine mode are always present. It will be interesting in fu- 
ture work to determine the extent to which the two BEC 
'mirrors' under consideration here can develop quantum 
correlations between the two light fields, leading e.g. to 
(a 2 ai) ^ 0, with similar behavior for the sine and cosine 
side modes. As such, this system may provide a sensitive 
test-bed to generate and study the onset of quantum cor- 
relations in this optically driven Bose condensate. These 
intriguing questions will be the topic of subsequent work. 
In the context of a BEC driven inside a Fabry-Perot res- 
onator, work has been done to treat quantum noise and 
correlations, see for instance [liUllES- Concentrating 
for now on the classical properties of the system, we turn 
in the next section to a discussion of its steady-state prop- 
erties. 



III. STEADY-STATE SOLUTIONS AND 
OPTICAL BISTABILITY 

Setting the time derivatives to zero in Eqs. (|5])- (fTTj) 
yields after some algebra a fifth-order polynomial equa- 
tion for |«i| 2 . This equation is quite cumbersome and 
is not presented here. We proceed by numerically deter- 
mining |ai| 2 , retaining only the physically relevant real 
positive roots of the fifth-order equation. From these val- 
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Figure 2: Mean intracavity intensity(unitless) of mode 1 or 2 
as a function of the detuning A for L — lOOnm, \ p = 780nm, 
K = 2tt x 1.3MHz, u) r = 2tt x 3.8kHz, [7 = 2tt x 3.1kHz, N = 
9000, 771=772 = 0.54k, and 7 = 0.001k. The dashed segments 
are unstable. The plot for modes 1 and 2 are identical, except 
that when mode 1 is on the CFD branch then mode 2 is on 
the CED branch, and conversely, see text. 



ues one readily finds |q!2 | 2 , X c , X s , ct\ and a 2 from 
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where 



and 



d = 8w r W 2 /(16w 2 + 7 2 ) 



C 2 
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The stability of the steady state solutions is determined 
by the standard Routh-Hurwitz criterion (38|. 

Figures [5] - [5] show typical results obtained for param- 
eters similar to those used in previous experimental Q 
and theoretical 0] work. For these specific parameters 
we obtain only three real roots for |ai| 2 . 

Consider first the case of real and equal pump ampli- 
tudes, 771 = r/2- Figure [2] shows the mean intracavity 
intensity of the counterpropagating optical fields 1 and 2 
as a function of the effective detuning A. Starting from 
large negative detunings, the intracavity field intensities 
increase as the magnitude of the detuning is decreased, 
as expected, with both field intensities being equal. At 
point A the system makes a discontinuous jump to point 

B. A similar behavior has been observed e.g. in the case 
of a BEC trapped in a Fabry-Perot cavity [J]. However, 
a new feature of the ring cavity system appears at point 

C. Here, the solution described by the thick curve be- 
comes unstable and the system undergoes a spontaneous 




Figure 3: Phase of the intracavity fields as a function of the 
detuning A for the same parameters as Fig. 2. The plots for 
modes 1 and 2 are identical, except that when mode 1 is on 
the upper semi- loop then mode 2 is on the lower semi-loop, 
and conversely, see text. The dashed segments are unstable. 
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Figure 4: Position X c (unitless) of the cosine mode as a func- 
tion of the detuning A. (a) Equal pumping amplitudes ni — 772; 
(b) slightly imbalanced pumping 771 = 0.53k; 772 = 0.54k; 
other parameters as those in Fig. [2] The dashed segments are 
unstable. 



symmetry breaking, with one of the intracavity intensi- 
ties increasing and following the CFD curve, while the 
other decreases and follows the CED curve. Which of the 
two counterpropagating fields will follow which branch is 
completely random and determined by classical (or quan- 
tum) noise, as illustrated in the following section. This 
behavior is reminiscent of that found for instance in sym- 
metrically pumped nonlinear interferometers [3ol [39| . 
Figure [U^a) shows the steady-state position X c of the 




Figure 5: Position X s (unitless) of the sine mode as a function 
of the detuning A for equal pumping amplitudes. Parameters 
as in Fig[2] The dashed segments are unstable. 
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Figure 7: Mean intracavity intensity(unitless) of mode 2 as 
a function of detuning A. Same parameters as those in Fig[6] 
The dashed segments are unstable. 




Figure 6: Mean intracavity intensity(unitless) of mode 1 as a 
function of detuning A. Parameters as those in Fig[2] except 
that r/i — 0.54k; and 772 = 0.55k. The dashed segments are 
unstable. 



cosine mode for the parameters of Fig. [5] Clearly, X c 
is also multistable, although this would appear not to 
be the case from the figure. This is because the branch 
CGD is actually degenerate, with the field intensities cor- 
responding to both branches CED and CFD in Fig. [2] 
yielding the same value of X c for a given detuning A. 
This might appear surprising, considering that these in- 
tensities are vastly different. The point is that the inten- 
sity only gives part of the story; we also need to consider 
the phase of the intracavity fields, see Fig. 3, which are 
just as important in determining the steady-state value 
of X c and X s , as apparent from Eqs. fl4")) and (JTSJ). More 
precisely, it is the interference between the forward and 
backward field amplitudes that drives X c and X s . To 
the left of the spontaneous symmetry breaking point C 
in Fig. [2] the phases of the two counterpropagating fields 
are equal, hence a^ai — oliol\ — 0, and the sine side 
mode is in a dark state, cf. Fig. 5. The situation is the 
same to the right of point D in Fig. 2. For detunings 
in the spontaneous symmetry breaking region, the upper 
intensity branch corresponds to the lower phase branch. 
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Figure 8: Position X s (unitless) of the sine mode as a func- 
tion of detuning A. Same parameters as those in Fig(6] The 
dashed segments are unstable. 

For a fixed detuning, these combinations of phase and 
amplitudes of the fields result in equal values of the in- 
terference term a^ai + ol-2.ol\ for the two branches, and 
hence the same value of X c . In this region we also have 
that a\ot\ — ol20l\ ^ 0, resulting in non-zero values for 
X s , see Fig.[SJ Note however that for that mode, the two 
branches lead to two values of X s exactly out of phase 
with each other, as expected from the form of Ea. (TT5|) . 

A small imbalance between the two (real) pump field 
amplitudes, 771 7^ 772 is sufficient to lift the degeneracy 
and remove the symmetry between the two counterprop- 
agating modes, see Fig. 4(b). This leads to a qualitatively 
different steady-state behavior of the system, including 
the appearance of isolated domains of solutions, or iso- 
las. A similar behavior has been previously predicted in 
other nonlinear optics contexts, including nonlinear ring 
resonators filled by a Kerr nonlinear medium [3~4l [3o| . 

Figure [5] shows the mean intracavity intensity of mode 
1 as a function of the detuning A for imbalanced pump- 
ing, and Fig.[7]the corresponding curves for mode 2. The 
effect of imbalanced pumping is to detach the loop pre- 
viously associated with spontaneous symmetry breaking 
from the "main" bistability branch. The isolated branch, 
the isola, is accessible via an adiabatic sweep of the de- 
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Figure 9: Time evolution of the intensity(unitless) of mode 1 
as the effective detuning is varied at a rate of 2ir x 0.3MHz 
from —4k to 4k. Due to random noise, the system evolves 
either along (a) the lower stable loop or (b) the upper stable 
loop, cf. Fig[2] The values given above the time axis are 
the instantaneous values of the detuning in units of k. Other 
parameters as in Fig[2] 



Figure 10: Time evolution of X s (unitless) as the effective 
detuning is swept linearly at a rate of 2n x 0.3MHz from —4k 
to 4k. Classical noise can result in the system evolving either 
along (a) the upper or (b) the lower stable loop, see Fig[5] 
The values given above the time axis are the instantaneous 
values of the detuning in units of k. Other parameters as in 
FigH 



tuning from large negative values. When reaching point 
C in Fig. 6 the main branch becomes unstable, and the 
field jumps toward point D on the isola. One could also 
reach the isola with a hard non-equilibrium excitation to 
a point within the basin of attraction of the steady-state 
branch RDE. That basin of attraction can be determined 
numerically by integrating the equations of motion (|8[)- 
(1111) backward in time from initial values arbitrarily close 
to that branch. 

Fig.BJb) shows X c as a function of detuning, illustrat- 
ing particularly clearly the lifting of the degeneracy asso- 
ciated with symmetric pumping and the onset of an isola, 
and Fig.[5]the corresponding X s . Further results not dis- 
cussed here show that when the difference between the 
pump field intensities are sufficiently large, both bistabil- 
ity and the isolas eventually disappear. 



IV. DYNAMICS 

We now turn to a discussion of selected aspects of the 
dynamics of the system as external parameters are var- 
ied in time. Here we concentrate on linear sweeps of the 
effective detuning from large negative values to large pos- 



itive values. Because the cavity field decay rate k is typ- 
ically orders of magnitude larger than the characteristic 
frequency of mechanical motion 4w r and the decay rate 
7 of the condensate, we can eliminate the optical field 
degrees of freedom adiabatically, and are left with a pair 
of coupled nonlinear second-order differential equations 
for X c and X s that can easily be solved numerically. The 
initial conditions for the oscillator variables are taken to 
be X c (0) = X s (0) = X c (0) = X s (0) = in all cases. 
We also add a small amount of classical noise 6N(t) to 
simulate atom number fluctuations, due e.g. to collisions 
or to thermal effects in the condensate. This noise term 
has the main effect of helping drive the system away from 
unstable branches in a reasonable time. In the simula- 
tions presented in this section 5N(t) is assumed to have 
a normal distribution with zero mean and the standard 
deviation equal to 5% of the total mean atom number. 

We consider first the case of real pumps of equal ampli- 
tudes, 771 =772, and sweep the effective detuning linearly 
at a rate of 2ir x 0.3MHz, a rate that is in most cases slow 
compared to the characteristic time scale of the atomic 
motion, in order to guarantee adiabaticity, but still fast 
compared to the atom loss rate [4j . Figure [5] shows the 
time evolution of the intensity of mode 1 as the detuning 
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Figure 11: Time evolution of the intensity(unitless) of mode 
1 under a linear sweep of the effective detuning at a rate of 
2tv x 0.3MHz, (a) for a sweep from —4k to 4k; (b) for a sweep 
from 4k to —4k, cf. Fig[6] The values given above the time 
axis are the instantaneous values of the detuning in units of 
k. Other parameters as in Fig[6] 



A is varied from —4k to 4k. Except for rapid transients 
that occur for a narrow range of detunings the inten- 
sity follows the steady-state values of Fig J3] These rapid 
transients correspond to system parameters such that the 
system is almost unstable, with a particularly slow relax- 
ation time to equilibrium. At the spontaneous symme- 
try breaking point C in Fig. 2, random noise determines 
whether the intensity follows the lower half loop solution 
of Fig. |H{a) or the upper half loop shown in Fig(5Jb). 
Figure HTJ1 shows the corresponding X s (t). 

We already mentioned that sweeping the effective de- 
tuning allows one to reach the isolas characteristic of 
imbalanced pumping. This is illustrated at point C of 
Fig. fTTTa). which shows the evolution of the intensity of 
mode 1 as the detuning is swept linearly from —4k to 
4k. Likewise, point P on Fig. [TTTb') shows that transi- 
tion when the detuning is swept down from 4k to —4k. 
For completeness, Fig. [T^] shows the corresponding time 
evolution of the position X s of the sine mode. 
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Figure 12: Time evolution of X B (unitless) as the effective 
detuning is swept linearly at a rate of 27r x 0.3MHz. (a) 
from —4k to 4k; (b) from 4k to —4k, cf. FigEJ The values 
given above the time axis are the instantaneous values of the 
detuning in units of k. Other parameters as in Fig[6] 



V. CONCLUSIONS 

In this paper we have presented a classical analysis 
of the fundamental optomechanical modes of a Bosc- 
Einstein condensate trapped inside a ring cavity. In con- 
trast to the situation of a Fabry-Perot cavity there are 
now two such modes, symmetric (cosine) and antisym- 
metric (sine) ones. These sine and cosine modes act as 
two coupled "condensate mirrors" of equal frequencies, 
their coupling resulting in a rich steady-state behavior, 
including for instance the appearance if isolas for an ap- 
propriate choice of parameters. 

One important feature of these modes is that together 
with the original condensate they form a V-system, with 
the upper levels - the sine and cosine modes - driven by 
a two-photon process involving both counterpropagating 
light fields. One or the other of these modes can then 
become a dark state as a result of destructive interfer- 
ences between the two counterpropagating fields. This is 
of course a classical effect: in the case of quantized opti- 
cal and matter- wave fields quantum fluctuations will nor- 
mally prevent these modes from becoming perfectly dark. 
It follows that measuring correlation functions of the op- 
tical field provides a direct means to probe the quantum 
properties of the matter-wave side modes. Future work 
will examine these aspects of ring-cavity optomechanics 
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in detail, both when quantum fluctuations can be treated 
as small fluctuations about a classical mean, and in the 
limit of very feeble fields where this approach is no longer 
appropriate. We will also revisit the Fabry-Perot case to 
include both the standing-wave mode that is normally 
included in the analysis, as well as the second frequency- 
degenerate but out-of-phase optical mode that needs to 
also be included and can become significantly excited for 
low mirror reflectivities. Finally, we will investigate in 
detail the damping of the matter-wave side modes due to 
the presence of a trap as well as the damping due to col- 
lisions. This will provide us with an effective Q-factor for 



the condensate mirrors and permit us to study in detail 
their heating and cooling mechanisms. 
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